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ABSTRACT 

We present a novel model in which field lines shortening after localized, three- 
dimensional reconnection heat the plasma as they compress it. The shortening pro- 
gresses away from the reconnection site at the Alfven speed, releasing magnetic energy 
and generating parallel, compressive flows. These flows, which are highly supersonic 
when P <^ 1, collide in a pair of strong gas-dynamic shocks at which both the mass 
density and temperature are raised. Reconnecting field lines initially differing by more 
that 100° can produce a concentrated knot of plasma hotter that 20 MK at the loop's 
apex, consistent with observations. In spite of these high temperatures, the shocks 
convert less than 10% of the liberated magnetic energy into heat — the rest remains as 
kinetic energy of bulk motion. These gas-dynamic shocks arise only when the recon- 
nection is impulsive and localized in all three dimensions; they are distinct from the 
slow magnetosonic shocks of the Petschek steady-state reconnection model. 

Subject headings: MHD — shock waves — Sun: flares 



1. Introduction 

Magnetic reco nnection has long been p roposed as a mechanism for heating coronal plasma. 



In one early model (|Kopp & Pneumanlll976l see Fig. [T]), reconnection occurs between open field 
lines separated by a vertical current sheet (red line) creating new closed field lines (post-flare 
loop s, grey). Closing these fie ld lines stops the solar wind upflow in a gas-dynar nic shock (GPS) 



that iKopp & PneumanI (|1976|) estimated would raise the temperature by 80%. ICargill & Priest 
(|1982l) found the direct magnetic energy conversion by reconnection to be a far more effective 
source of heating in this same model. It could raise the temperature of post-flare loops by up to 
a factor of three to 6 MK. Even this higher value is, however, insufficient to explain the 15-20 
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side view 

Fig. 1. — The classical Kopp-Pneuman model and its modifications. The side view (left) shows 
the basic two-dimensional geometry. Inset (a) on the upper left shows how steady Petschek recon- 
nection is configured within the current sheet. Inset (b) on the lower left shows how a transient 
reconnection produces a finite layer of reconnected flux (grey). The front view (right) shows how 
the model becomes three-dimensional when the reconnection is localized in third dimension (x). 
Instead of a current layer the reconnection produces a A-shaped flux tube of closed field (grey); a 
completely detached, V-shaped flux tube (not shown) is created at the same time. 
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MK temperatures o bserved at apices of post-flare loops in Soft X-ray and Fe XXIV EUV emission 
Jwarren et all 19991) . 



Recent theoretical investigations have revealed that electric fields large enough for fast mag 
netic reconnection can be self-consistently produced by a wide range of small-s cale mechanisms 



provided only that they a re localized within a segment of the current sheet (|Birn et alj 12001 



of lPetschekl (|1964[ see inset a) with slow magnetosonic shocks outside the non-ideal region (SMSs, 
blue lines) at which the magnetic field is deflected and weakened, thereby heating the plasma. 
When reconnection is localized in both space and time the S MSs close back together across a fi 



Biskamp & Schwarzll200lh . All such processes generate reconnection flows resembling the model 



nite l ayer of horizontal field, as shown in inset (b) of Fig. [Tl (|Semenov et al.lll983l : iBiemat et al 



19871). This field forms the top of a "hairpin" comprising all the flux closed since the onset of 
reconnection. 

As the hairpin flux sheet retracts its field lines become much shorter and releases substantial 
magnetic energy. This magnetic energy is converted al most entirely into kin etic energy rather 
than partly into heat as in Petschek's steady state model (iSemenov et al.ll 1998b . In a strictly two- 
dimensional version of the model the mass of the shortened tube accumulates in the tip of the the 
retracting hairpin. This "snowplowing" artifact is absent when reconnection is also localized in the 
third dimension (x) and there is a magnetic field component in that direction. With the horizontal 
field component (sometimes called a "guide field") the current sheet separates field with angle 
(tanC = BJB.^, see the "front view" of Fig. [B- 

Flux reconnected within a patch and ove r a finite interval forrns a A-shaped flux tube (grey ) 
similar to magneto spheric yZwx transfer events ( Russell & Elphic|[l978 : iLee et al.l[l993l : [ottoll 19951) . 
The tube is distinguished from the surrounding field (the flux layers) by the distinct connectivity 
given it through reconnection. Different field line geometry, such as the bend, produce dynamics 
in the tube entirely different from those in the surrounding flux layers, with which it has little 
subsequent interaction. 



Recently iLinton & Longcopd (120061) studied the relaxation of this post-reconnection flux tube 
using three-dimensional magnetohydrodynamic (MHD) simulation. They found that the perpen- 
dicular dynamics of the tube, as it mov es between t he flux layers, is well approximated by the 
equations for a thin magnetic flux tube (|Spruitill98l|) . although with P < 1. They presented an 
analytic solution in which the post- reconnection tube shortens at the Alfven speed, converting 
magnetic energy into kinetic energy. iLinton & Longcopd (120061) did not investigate the dynamics 
parallel to the tube or any associated thermal effects of shortening. 



This letter demonstrates that flux tube shortening is a powerful, inevitable mechanism for 
heating post-flare loops, which has not been previously investigated. The magnetic forces respon- 
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sible for shortening, also drive compressive parallel flows at the Alfven speed. At very low (3, these 
are high-Mach number flows whose collision naturally generates very strong shocks. The shocks 
are distinct from the SMSs of Petschek reconnection and are driven by reconnection-initiated, per- 
pendicular dynamics. This is in contrast to previous investigations of flux tube shocks wherein 



acoustic wave-steepening or pressure differences were considered as drivers (IHerbold et al.lll985 



Thomas & Montesinoslll991 ). 



2. Post-reconnection flux tube dynamics 

We begin by assuming that localized, transient, fast magnetic reconnection has occurred, by 
an unspecified physical mechanism, within an otherwise static current sheet. This will, as just 
discussed, leave a A-shaped flux tube, initially at rest. Due to its sharp bend it is out of equilibrium, 
and magnetic forces start it sliding downward between the magnetic layers separated by the current 
sheet. 

The tube's retracti on, unhindere d by the external flux layers, can be modeled using the thin 



flux tube equation s of iSpruiti (119811) and subsequent authors (see, for example, the review by 
Fisher et al.ll2000|) . While the high-/? flux tubes in those previous investigations are confined by 
the pressure of the unmagnetized convection zone, our post-reconnection tube has very low /? 
and is confined by the magnetic pressure of the flux layers outside the curr ent sheet. The tube 



is an isolated entity dist inguished from its surroundings by its connectivity (|Linton & Longcope 



2006l : lLinton et al.ll2008|) . We assume sufficient coUisionality to justify the use of MHD equations 



throughout. 

The tube is assumed thin enough to be described only by its axis. Internal properties such 
as the magnetic field strength, Bi, pressure, Pi, and mass density pi, are function only of axial 
position. The tube is also thin enough for fast magnetosonic waves to establish pressure balance 
across its diameter virtually instantaneously. This assumption constrains the internal properties to 
match those outside the flux tube: Bf /Sir + pi = B'^/Sn + pe, assumed to be uniform and constant. 

We will also assume the plasma P to be always small. The main force on a section of tube is 
therefore the magnetic tension due to curvature of the axis Bf{dt/di) /An, where i is arc-length 
and t = dx/diis the unit tangent vector. Since it is ultimately the Lorentz force, it is natural that 
this force is strictly perpendicular to the axis (i ■ di / d£ = 0). 

The pressure gradient, —idpi/di, is formally smaller, by a factor of than the magnetic 
tension. It is, however, the only force parallel to the axis, and is essential to arresting internal 
pile-up of mass. We therefore retain terms involving pressure to first order in /3. 
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Velocity evolution is governed by a momentum equation with only these two forces. Con- 
servation properties become apparent when arc-length is replaced with /i, the integrated mass per 
unit flux: d^ijdt = p j Bi. The value of fi never changes for a given fluid element. The resulting 
momentum equation 



dt 4:7rpi d£ \ Att 

includes the parallel pressure gradient after use of pressure balance, Bi = 3^, + 47r(pe — pi)/ B^, 
valid to first order in (3. The momentum per unit flux of any section of tube, 

P-^M = Jvdp , (2) 

changes only through forces (per unit flux) from the ends of the section, Bit/An, directed parallel 
to the axis. 



3. Shock relations in thin flux tubes 

Momentum conservation leads to a set of shock relations for thin flux tubes. Consider two 
straight sections with uniform properties (designated 1 and 2), separated by an abrupt change at 
coordinate po. The length-scale of this change is large compared to the tube radius but otherwise 
small enough that we hereafter call it a "discontinuity" and "comer"Q This feature moves through 
space at constant velocity u while its Lagrangian coordinate changes at constant rate fiQ. We 
assume that the properties of the straight sections do not change as this happens. Fluid velocities 
on either side of the discontinuity differ from u, but components perpendicular to the tangent vector 
match that of u: (v — u) x t = 0. 

The component of relative velocity parallel to the tangent vector, u ^ = t ■ (v — u), represents a 
flow across the discontinuity. The mass flux (per magnetic flux) through the discontinuity, /iq, must 
equal the mass flux across points on either side of it, since the comer moves without changing. For 
positions pi and 112, separated by fixed distances from po this means 

fi2 = - PivJBi =-piV jBi = fi^ = fiQ . (3) 



Analysis of the discontinuity's internal structure must be done outside the thin-tube approximation. Doing so will 
probably reveal its length-scale to be comparable to the tube's radius. This is analogous to non-ideal effects resolving 
internal structure of a hydrodynamic shock at scales comparable to the mean free path JGradlll95ll) . Following that 
analogy, the thin flux tube equations provide external conservation laws leading to the shock relations. 
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The flux tube between /i2 and yUi does not change so its momentum (per magnetic flux) is 
constant, 

P = l/ivl + Ii?.t/47r] = , (4) 

denoting [/] = /2 — /i- Since v = u + f ||t away from the discontinuity momentum constancy 
may be written 

(i?,/47r - p,vl/B,) |t] - Ip,/B, + p,vyB,\ (t) = , (5) 

where (/) = + /i)/2. 

The vectors [t] and (t) are orthogonal so eq. ([5]) represents two independent equations. As 
long as the bend is not a hairpin ((t) ^ 0) the (t) component of momentum conservation, to lowest 
order in (3, is equivalent to gas-dynamic momentum conservation: \pi + pjt>^] = 0. To the same 
order in /5, eq. ([3]) gives the gas-dynamic mass continuity equation: IpifJ = 0. Provided there 
is no heat flow from the external field or across pi or then the sum of kinetic and thermal 
energy within the tube section will also be conserved. This provides one additional, independent 
constraint on pi, pi and v^^ across the discontinuity. (Since Bi can be related directly to pi it can be 
eliminated from energy conservation.) 

The foregoing describes three relations between six quantities which do not include the mag- 
netic field' s strength or direction. The r elations are the traditional gas - dynamic Rankine-Hugoniot 



conditions ( Courant & Friedrichslll948 ). but for flows inside a flux tube ( Ferriz-Mas & Moreno-Insertis 



1987|) . The other component of eq. dS]), to lowest order in P, 

(l-(4Wi?eX)lKl' = ■ (6) 
constitutes one more relation, which does involve the field direction. 

From values on one side the Rankine- Hugoniot relations may be s atisfied in two different 



ways by values on the other (see for example ICourant & Friedrichslll948 ). They may be satisfied 



non-trivially by a unique set of dijferent values, or they may be satisfied trivially by the same set 
of values (i.e. |pj] = It'll ] = IpJ = 0). In the non-trivial case, the three different quantities satisfy 
three independent constraints and cannot be forced, in general, to satisfy a fourth. The leading 
factor of eq. however, constitutes a fourth independent constraint so it will not in general 
vanish; it is therefore necessary that [t] = (there is no bend). This means that a thin, low-/5 flux 
tube can support discontinuities in internal quantities only at a GDS within a straight section of 
tube. 

If, on the other hand, there is a bend in the flux tube (|t] ^ 0), the only way to satisfy eq. © is 
for the pre-factor to vanish. Since the non-trivial solution of all three Rankine-Hugoniot conditions 
would over-determine the system, they must be satisfied trivially, without discontinuity. In other 
words Pi, and v are continuous at pq and satisfy the relation v = B^/ y^Anpi = va the Alfven 
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speed. This is similar to an intermediate shock (IPriest & ForbesI 120001) . but includes the influence 



of fast magnetosonic waves assumed to be maintaining pressure balance across the tube. 



4. Shocks in the retracting flux tube 

The A-shaped bend in the post-reconnection flux tube shown in Fig. [T] will immediately de- 
compose into four different shocks of the kinds described above. Two intermediate shocks (bends, 
B) propagate along the field lines, forming a straig ht horizontal section between them (see Fig. [2l), 



as previously found by Linton & Longcope ( 2006|) . Two GDSs propagate away from the center. 



at ±Vs, along the horizontal section. This symmetric arrangement divides the tube into sections 
labeled, 3, 2 and 1, outward from the center. 

Section 1 consists of the flux tube at rest and in its initial state: vi = 0, pi = pe and pi = Pe- 
The initial field lines are inclined at angle (, so ti = x cos C — z sin (, while t2 = x in the horizontal 
section. The bend propagates along the initial field at the Alfven speed: u = f ^.e^i- The relative 
parallel flow is continuous across the bend t>|| i = f y 2 = — ^A,e- These facts can be combined into 
the post-bend fluid velocity 

V2 = u + t;||t2 = -2t;^^eSin^(C/2)x-t;^,eSinCz , (7) 

directed along the bisector of the bend (the mean direction of the curvature force). 

Internal quantities are continuous across the bends, so P2 = Pi = Pe and P2 = Pi = Pe- From 
a reference frame moving downward at — zf^ ^sin^, sections 2 and 3 appear to form a classic 
shock tube, with inflow at 

M2 = ^ = ./^sin^(C/2) , (8) 

where 7 = 5/3 is the ratio of specific heats. Since /3e <^ 1, this Alfvenic inflow can have extremely 
high Mach number when field lines of significantly different orientation reconnect. 

This supersonic inflow is brought to rest, Vx,3 = 0, by a GDS moving outward at Vg. This 
stopping shock, equiv alent to a piston moving in to stationary fluid at speed \va:^2\, is a classic 



problem (see §69 from lCourant & Friedrichsll 19481 for example), for which the solution is 



Vs/K2\ = ^M2-2 + (7 + 1)2/16 - K3 - 7) • (9) 

The density ratio, P3/P2 = 1 + \vx,2\/vs, following from mass conservation, approaches the well- 
known limit ps/pe = (7 + l)/(7 — 1) = 4, at large Mach numbers. 

The post-shock pressure, = /5e[l + 7^|(1 + Vs/\vx,2\)], also follows from the shock 
relations, and even in the limit of vanishing pre-shock pressure (Pe ^ 1)^ it can be significant: 
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GDS GDS 

Fig. 2. — The flux tube retracting from the A-shaped initial condition shown in Fig.[il This initial 
configuration is shown by thick dashed lines, while its current state is shown in grey. The past 
trajectory of the fluid, shown by solid lines, is initially focused inward. The gas-dynamic shocks 
(GDS) redirect the flow vertically; their past positions are shown by a thin dashed line. Circled 
numbers designate the states, 1, 2 and 3, separated by these shocks. 
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/?3 — 4(7 + 1) sin'^(C/2). The low-/5 assumption thereby imposes a limit on the reconnection 
angle; only for C < 36° is /^g < 0.1 and /Jg > 1 when C > 67°. 

The post-shock temperature is more conveniently expressed with respect to = {nip/ ^ e/2, 
than to the initial temperature. For example, 15 G field immersed in rig = 3 x 10* cm~'^ plasma has 
a characteristic temperature = 2 x 10* K. The plasma beta is Pe = Te/T^, so a 2 MK coronal 
plasma will have Pe = 10^^. Figure [3^-b shows the post-shock temperature, T3 = T^jS^Pe/ Pz, 
over a range of Pe-C, parameter space. For /3e = 0.01 and /^s = 1/3 (the edge of the light grey area, 
C ~ 50°) the post-shock plasma will be T3 = O.IT^^ = 20 MK. 

The total energy release initiated by reconnection, = 2v^ ^Lpe sin^(^/2), is the initial 
kinetic energy of all mass within the L = 2^;^ gt of tube affected by retraction. This is greater than 
the energy decrease resulting from shortening the tube by AL = 2L sin^((^/2), due to additional 
work done by the background magnetic field expanding into the vacated volume. In our idealized 
model this work exactly doubles the energy, but in more realistic scenarios the factor may be 
somewhat different. 

The moving mass is deflected downward by the GDS, converting a portion of its kinetic 
energy into thermal energy. The kinetic energy (per area) converted to thermal energy, p^vl 2Vst, 
constitutes a fraction 

|| = JS!|!^ = 2(l + ^)sin^C/2). (10) 

of the total (see Fig. [3]:). Cases of vanishing initial pressure (/3e <C 1) will have a fraction (7 + 
1) sin'^(^/2) of the released energy thermalized. Compression work done on larger initial pressure 
will raise this fraction slightly, but all cases with /?3 < 1/3 thermalize less than 10% of the released 
energy. 



5. Summary 

The foregoing illustrates, through a simplified analytic model, a process we believe must 
be common in the flaring corona. Localized reconnection does not directly dissipate magnetic 
energy, but rather initiates its release through subsequent shortening of field lines. This shortening 
propagates from the reconnection site at the Alfven speed, converting energy from magnetic into 
kinetic form. The shortening flux tubes compress plasma within them, raising its temperature as 
they do so. Due to the supersonic (Alfvenic) flows this compression occurs at strong shocks, which 
raise the temperature far beyond that from adiabatic compression. 

Both the GDS in our model and the SMS in Petshcek's model result, ultimately, from shorten- 
ing or weakening of magnetic field lines. The two-dimensional Petschek model permits shortening 
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only perpendicular to the symmetry direction, accompanied by weakening. The shock is therefore 
a SMS whose normal is mostly perpendicular to the sheet. After three-dimensional patchy recon- 
nection, on the other hand, field lines also shorten in the erstwhile symmetry direction which is 
the orientation of the shock normal (±x). The GDSs are disconnected from the diffusion region, 
and are instead features of the ideal relaxation following reconnection. This new scenario requires 
two-step energy conversion, from magnetic to kinetic to thermal energy, in contrast to the SMSs 
which thermalize magnetic energy directly where they weaken the field. 

This work was supported by NASA and the NSF. 

REFERENCES 

Biemat, H. K., Heyn, M. R, & Semenov, V. S. 1987, JGR, 92, 3392 
Bim, J., et al. 2001, JGR, 106, 3715 
Biskamp, D., & Schwarz, E. 2001, Phys. Plasmas, 8, 4729 
Cargill, R J., & Priest, E. R. 1982, Solar Phys., 76, 357 

Courant, R., & Friedrichs, K. O. 1948, Supersonic Hows and Shock Waves (New York: Inter- 
science Publishers) 

Ferriz-Mas, A., & Moreno-Insertis, F. 1987, A&A, 179, 268 

Fisher, G. H., Fan, Y., Longcope, D. W., Linton, M. G., & Abbett, W. P 2000, Phys. Plasmas, 7, 
2173 

Grad, H. 1951, Comm. in Pure and App. Math., 5, 257 

Herbold, G., Ulmschneider, P, Spruit, H. C, & Rosnr, R. 1985, A&A, 145, 157 
Kopp, R. A., & Pneuman, G. W. 1976, Solar Phys., 50, 85 
Lee, L. C, Ma, Z. W., Fu, Z. F, & Otto, A. 1993, JGR, 98, 3943 

Linton, M. G., DeVore, C. R., & Longcope, D. W. 2008, Earth, Planets and Space (in press) 
Linton, M. G., & Longcope, D. W. 2006, ApJ, 642, 1177 
Otto, A. 1995, JGR, 100, 11863 



-11- 



Petschek, H. E. 1964, in AAS-NASA Symposium on the Physics of Solar Flares, ed. W. N. Hess 
(Washington, DC: NASA), 425 

Priest, E. R., & Forbes, T. G. 2000, Magnetic Reconnection. MHD theory and applications (Cam- 
bridge University Press) 

Russell, C. T., & Elphic, R. C. 1978, Space Sci. Rev., 22, 681 

Semenov, V. S., Heyn, M. F, & Kubyshkin, I. V. 1983, Sov. Astron., 27, 600 

Semenov, V. S., Volkonskaya, N. N., & Biemat, H. K. 1998, Phys. Plasmas, 5, 3242 

Spruit, H. C. 1981, A&A, 98, 155 

Thomas, J. H., & Montesinos, B. 1991, ApJ, 375, 404 

Warren, H. P., Bookbinder, J. A., Forbes, T. G., Golub, L., Hudson, H. S., Reeves, K., & Warshall, 
A. 1999, ApJ, 527, L121 



This preprint was prepared with the AAS L^TeX macros v5.2. 



-12- 




Fig. 3. — The post-shock temperature and thermal energy across the (-(^e parameter space. The 
top panel (a) shows contours of T3 whose levels can be determined from intersection with the left 
axis (C = 0) since T3 = there. Solid contours show T3 = 0.01, 0.1 and 1 (all in units of T^). 
Light and dark grey regions show 0.333 < P3 < 1 and > 1 respectively. The middle panel (b) 
shows T3 versus C for the values = /3e = 0.01 (squares) and 0.1 (triangles). The bottom dashed 
curve is limit for the case Pe — 0. The bottom panel (c) shows, for the same values of = /3e, the 
fraction of released energy thermalized by the GDS. Symbols show the same point in each panel. 



